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Qh! Abstract 

r-| ! We re-examine like-sign dilepton signals at the Fermilab Tevatron assum- 

ing that the excess high-Q 2 events recently seen at HERA are due to scalar 
resonances such as squarks of i?-parity-violating supersymmetry. For gluinos 
in the mass range of 200-350 GeV, the like-sign dilepton signal can help to 
make the crucial distinction between the most favoured squark explanation 
and other proposed solutions. 
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The particle physics community has recently been intrigued by reports of excess frigh- 
ts 2 events in both the HI and the ZEUS detectors at HERA. Concentrated near a 
parton- momentum fraction of x ~ 0.5, these events appear to point to some new effects 0- 
n[| beyond those expected in the Standard Model (SM). Although the contrary opinion 
has also been expressed jR], the jury is still out. Even assuming — tentatively, in view 
of the low statistics — that new physics has indeed been found, opinions differ as to what 
these effects could be. Various suggestions have been made, including compositeness @, 
contact interactions [RPJ, and leptoquark or squark pj-fij resonances around 200 

GeV. Of these, the last - - namely squark resonances - - constitutes one of the more 
exciting possibilities, since these events could then be the first experimental signals of 
supersymmetry (SUSY) — a desirable option for various theoretical reasons ITU 



The possibility of a squark resonance at HERA was first considered by Hewett [14] 



and subsequently by others fl5| . This obviously requires a violation of i?-parity |16| in 
the form of lepton-number- violating A' operators []H]]. The relevant Lagrangian is given 
by 



C\> — —Xij k 



ViLdkRdjL + dj L d kR v iL + dl^fiLYdjL 

(1) 



-ZihdkRUjL - u jL d kR e iL - dl R (e iL ) c u jL 



+ h.c. 



where the effects of quark mixing have been neglected. In a previous work Q, the present 
authors pointed out that, if the excess events at HERA are caused by a single dominant 
A' coupling, then the relevant ones are A' 121 (with production of a cl squark), or A' 131 (with 
production of a t L squark). A third possibility has also been suggested |5],|] — involving 
A / 132 , with a ti resonance. A value of A' m in the interval 0.03-0.26 could give rise to the 
appropriate signal for the first two options. The third would require A' 132 close to its ex- 
perimental upper bound of 0.3. Some other operators can contribute, though marginally, 
for the highest experimentally allowed values of the relevant coupling constant j5]. 



A crucial feature of the i?-parity- violating (^t p ) signal is the fact, emphasized by all 
workers in this field, that the squark can also decay through i?-parity-conserving channels 
to charginos and neutralinos, which could give rise to distinctive and unambiguous signals 
at HERA |18|]. Observation of such a signal will be an immediate confirmation of the 
SUSY hypothesis. We understand that an intensive study of this nature is indeed under 
way [jl9| . However, non-observation of such signals cannot rule out supersymmetry — it 
will merely tell us that the direct _R-parity-violating decay of the squark to e + + d has 
to be the dominant decay channel, while gaugino decay channels are suppressed by large 
gaugino masses or small couplings. In this case, of course, the decays of the squark will 
be identical with those of a leptoquark with the same quantum numbers. 



An interesting feature of the -R-parity-violating SUSY solution is that low values of 
the A' coupling in question require large branching ratios to the .R-parity-violating decay 
channel un — > e + + d. While this may seem paradoxical at first, we need to recall that the 
number of excess events seen at HERA essentially scales as X' 2 /3, where f3 is the branching 
ratio of un to e + + d. Thus, small values of A' require a large branching ratio and vice 
versa. Such a scenario is far from unnatural — it can be achieved without difficulty by 
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requiring the charginos and neutralinos to be either heavy or Higgsino-like, which merely 
corresponds to a restricted range in the parameter space of the minimal supersymmetric 
Standard Model (MSSM). 



In view of the above, it is now amusing to ask the following question. If, with more 
statistics from the 1997 run, the high-Q 2 anomaly does survive as a genuine effect and, 
moreover, as one compatible with a particle resonance, but no further signals of SUSY are 
seen at HERA, then the resources of that experiment will have been exhausted insofar 
as telling the nature of the new physics is concerned. It could be a leptoquark or a 
squark with a small ( 0.1) coupling. In this case, can any other running (or projected) 
experiment distinguish the supersymmetric option from the others? 



The only other running high-energy facility with sufficient energy to produce the 
'particle' under consideration is the Fermilab Tevatron, running at a centre-of-mass energy 
of 1.8 TeV. Squarks or leptoquarks can be pair-produced at this machine and their decays 
would give rise to distinct dilepton + multijet signals. For a squark mass of 200 GeV, we 
estimate the pair production cross section to be around 0.19 pb, which is consistent with 
Ref. |20fl. Searches for such processes by the DO Collaboration have already put a bound 



of 175 GeV (with j3 ~ 1) on the mass of the resonance at 95% confidence level ||21|| . This 



bound is based on a data sample of 117.7 pb -1 and may be expected to increase to cover 
the mass range of interest (~ 180-220 GeV) as more data accumulates. If these searches 
fail to show up a resonance, it will be a sign of some serious glitch(es) in our understanding 
of the HERA anomaly, since most of the new physics solutions — including SUSY — do 
predict corresponding signals at the Tevatron. On the other hand, if the dilepton + jets 
signal is confirmed, it will inevitably lead us to the question of distinguishability between 
various theoretical models. 



This report addresses the crucial question of distinguishability for the SUSY solution 
in the context of Fermilab Tevatron searches. An interesting feature of SUSY is that the 
complete theory requires many more (light ?) particles than the squark. One of these is 
the gluino — the fermionic superpartner of the gluon — which is a strongly interacting 
Majorana fermion. As a result, relatively large cross sections from QCD production of 
gluinos are possible and processes involving them can violate fermion number, leading to 
rather spectacular like-sign dilepton signatures. These have already been studied in some 
detail in the literature p2|-|2~4}| and our main purpose is to re-examine their utility in the 
context of the squarks that can explain the HERA events. The importance of this study 
lies in the fact that such signals will be present only if the HERA excess is caused by a 
squark resonance. 



We thus need to consider the pair-production of gluinos at the Fermilab Tevatron, 
followed by the decay of each gluino into a quark and a squark, where the squark is 
the same one that explains the HERA anomaly. Because of the Majorana nature of the 
gluinos, the squarks coming from their decay have equal probability of forming a like-sign 
or an unlike-sign pair. Subsequent decays of these squarks to e + + d (e~ + d) would lead 
to equal numbers of like-sign and unlike-sign dilepton pairs. This decay chain differs from 
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those studied in Refs. [^2],f23|] in that the squark now has no decays through gaugino-led 
channels, the latter being the exclusive decay modes studied therein. 



Once produced, the gluino will decay into all possible quark-squark pairs that are 
kinematically allowed. Concentrating on the (A' 121 ,cl) scenario — the most promising one 
in the context of the HERA events — the decay channels of interest are 

g ^ c + cl — > c + e + + d (2a) 

and 

g — > c + c* L — > c + e~ + d , (2b) 

which leads to a final state with a dilepton pair and jets. The combinations e + e + , e + e~, 
e~e~ appear with probability ratio 1:2:1. In other words, if Cl were the only light squark, 
the decay of a gluino pair would lead to a like-sign dilepton final state in exactly half the 
cases. 



Now, the branching ratio of the gluino to e + + X is given by 

B ^ r( ^ e++X)= ^(^)W , (3) 

where Xi = m^jM?, X{ = m^./M- and the sum over i encompasses all squarks for which 
the decay g — ► is kinematically possible. In Eq.(|^), (3i is the branching fraction of 



the relevant squark to e + + X, while A (a, b, c) = J (a — b — c) 2 — 46c. From our previous 
discussion, we see that, for A ^ 0.1, 

As long as A' 121 is the only non-zero Ifl p coupling, the cascade deacys of the other squarks 
will result in final states with or (anti-)neutrinos. Such event topologies have been 
considered in Refs. [22], 23]. Rather than repeat their analyses, which are quite compre- 



hensive for these decay chains, we shall make the conservative assumption that — 
for these squarks. With this simplification, B in Eq. (§) is determined solely by the gluino 
and the squark masses. We find that for A' 121 ~ 0.1, B is 0.54 or less. 



As the masses of the squarks are unknown parameters in the MSSM, we have to 
make some assumptions for the sake of simplicity. An explanation of the HERA anomaly 
requires only two squarks to be light. These are the cl (obviously) and the §l whose mass 
is given by the relation 

m? L = m 2 - L — m 2 c +m 2 s — cos 2f3 . (4) 

Numerically, for m c - L ~ 200 GeV, we find that varies in the range 200-215 GeV as 
tan/3 varies from 1 to 50. The other squarks can be as heavy (or as light) as we please. 
We make the ad hoc assumption that 

m um = m d m = m u La = m q > ( 5 ) 
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where i = 1, 2, 3 and a = 1,3. The masses of di and hi are given by relations analogous to 
that in Eq.(|j). This is not perhaps the optimal spectrum for the signal in question, nor is 
it the one most favoured by model-builders, but it is certainly the simplest one. We have 
checked that the qualitative features of our results do not change with different Ansatze 



for the squark masses, although the detailed numerics are somewhat affected [25 
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Figure 1: Contours for (3c L = 0.9 in the M^-\x plane for tan/5 = 2 and various values of 
the ratio M\jMi. Dashed, solid and dot-dashed lines correspond to M1/M2 = 2,0.5,0.25 
respectively. The dependence on tan (3 is not very significant. 

Since we want the lfl p channel to dominate the squark decay width, the gluino must 
be heavier than the squark. And so must be the charginos and neutralinos, unless of 
course the lighter ones are Higgsino-dominated. Within the framework of gaugino mass 
unification, and the LEP constraints, this implies a large value for the gaugino mass 
parameter Mi. Within the same framework, this would further imply very large gluino 
masses (generally ^ 500 GeV) and hence extremely tiny rates at the Tevatron. Note, 
however, that gaugino mass unification is not very well explored within the context of im- 
parity violation. While Grand Unification can perhaps be made consistent with Jji v f26| , 
such models typically need extra multiplets, all consequences of which are not very well 
studied. Hence, we shall not assume gaugino-mass unification. The squark mass spectrum 
does not assume scalar mass unification either. We thus adopt a purely phenomenological 
approach to the SUSY mass spectrum. 



Once we give up gaugino mass unification as a hypothesis, the soft supersymmetry- 
breaking parameters Mi, M2, M$ of the MSSM are no longer related. It is then possible 
to obtain large branching ratios for cl — > e + + d for a fairly large range of the parameter 
space. This is illustrated in the Mi~\x plane for a fixed value tan j3 = 2 in Fig. fig: branching. 
We have chosen two values of the coupling (a) A' m = 0.05 and (b) A' 121 = 0.1, both of 
which require a large branching ratio in order to explain the HERA events (see Fig. 2 
of Ref. [[J]). The region above each contour corresponds to j3c L > 0.9. Dashed, solid 
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and dot-dashed lines correspond to the ratios Mi /Ma = 2, 0.5, 0.25 respectively, of which 
the second approximates closely to the ratio M\jMi = §(<7i/<72) 2 obtained by assuming 
gaugino mass unification purely in the electroweak sector of the MSSM. It is also worth 
mentioning that for this figure we have included LEP constraints on the 1/lp scenario in 
question — these come only from the total width of the Z, to which the excess SUSY 
contribution is ]2"7fl : 



£ r(Z - xtxj) + E m - X°X°) S 23.1 MeV . (6) 

i,j=l i,j=l 

It is at once apparent that somewhat large values of M 2 (and Mi) are demanded. Low 
values of are disallowed, partly by the LEP constraint and partly because they cor- 
respond to small branching ratios. Finally we should note that varying tan/3 does not 
change the contours much. Accordingly, we have fixed tan/? = 2. Our plots are in good 
agreement with those of Ref . . 



The like-sign dilepton signal has been estimated assuming (3g L = 1 and using a parton- 
level Monte Carlo event generator. Our estimates for gluino pair-production agree well 
with those of Ref. [pOR . For our simulation, we use the CTEQ3M structure functions p8| , 



which were calculated using the package PDFLIB We assume a detector rapidity 

coverage of 3 for both jets and leptons : 

m\,\vm<z- (7) 



All putative jets that lie within a cone of AR = \J Ar] 2 + A</> 2 < 0.7 of each other are 
merged (the momenta added vectorially) to form a single jet. Here Ar) is the difference of 
their rapidities and A<j) denotes their azimuthal separation. The major SM backgrounds 
to the like-sign dilepton signal emanate from ti and WZ production |22|, p3f . These can 



be severely suppressed by using suitable kinematical cuts. We demand that the leptons 
have sufficient transverse momentum and be relatively isolated : 

p T (£) > 15 GeV, E* c (£) < p T (£)/A , (8) 

where E^(£) denotes the total hadronic energy within a cone of AR = 0.4 around the 
lepton £ . With the imposition of these cuts, the SM background due to ti production 
was estimated to be as low as 0.3 fb, while 2.1 fb come from WZ production P2"fl . Clearly, 
the background level is very low and thus the discovery limit is primarily governed by the 
signal strength. Our results are consistent with those of Ref. | ]23|] , which were checked 
against an ISAJET-based calculation with reliable cuts and efficiencies. 



In Fig. we present contours of cross-section for the like-sign dilepton signal in the 
plane of rrig, the common squark mass, and M g , the gluino mass. We have fixed the mass of 
the cl at 200 GeV and tan j3 = 2 (which fixes the mass of sl ~ 210 GeV). We have checked 
that a variation of the former between 180 and 220 GeV and the latter between 1 and 50 
makes little visual impact on the figure, although the detailed numbers do change slightly. 
It should be noted, though, that the lower edge of the curves represent conservative 
estimates as we have neglected the cascade decays of the squarks into the dilepton channel. 
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Figure 2: Contours for the like-sign dilepton cross-section at the Tevatron in the m- g ^mq 
plane. The dependence on tan (3 is minimal. Changes in m 5L affect only the lower edge 
of the curve (see text). 



Inclusion of these is expected to make the low-m^ region more accessible [22,23]. Dotted 



lines in the figure indicate the kinematic reaches of the Tevatron with 0.1 fb _1 and 1 fb -1 
of integrated luminosity. In view of the low background of ~ 2.4 fb, we show contours 
of 5 fb and 50 fb, since these imply 5 events at the Tevatron with, respectively, 0.1 fb" 1 
of data (which are currently available) and with 1 fb -1 of data (which should become 
available with the Main Injector run). The solid contours may, therefore, be taken as 
discovery limits. The dot-dashed line corresponds to a signal size of 200 fb, which is 
perhaps rather optimistic. It is immediately obvious that if Nature favours us with a 
somewhat low-lying gluino state (< 350 GeV), then it should be possible to observe like- 
sign dilepton signals at the Tevatron, unless, indeed, there also happens to be a bunch of 
low-lying squark states around 200 GeV apart from the (cl,§l) pair which have $ = 
(thereby reducing B). On the other hand, if these signals are not seen, it should be 
possible to rule out a substantial region in the squark-gluino mass plane, which has not 
been hitherto accessible to Tevatron searches (assuming, of course, that the HERA events 
are susceptible to a SUSY explanation with a cl resonance). 



We now comment on the equally interesting possibility that the HERA events are due 
to a stop t L resonance through a A' 131 (or A' 121 ) coupling. In this case, a strong like-sign 
dilepton signal requires the gluino to decay to a top-stop pair which is only possible when 
the gluino mass is about 370 GeV or more. This mass range is at the kinematic limit of 
the Tevatron and, as Fig. |2] makes clear, yields a very small signal. Thus, we conclude 
that the top-stop decay mode of the gluino will not yield an observable like-sign dilepton 
signal through the direct Ifi p decay mode. However, the story does not end here, because 
the presence of an ~ 200 GeV ti implies the presence of a light hi. Recalling that the 
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masses are related by 



= m\ L — m\ + m\ — m\r cos 2(3 , 

it is easy to check that a 200 GeV ti implies a &l in the mass range 105-132 GeV. (The 
lower part of this range is ruled out in the limit of vanishing left-right stop mixing by 
constraints from the p parameter, although in the presence of mixing this constraint can 
be relaxed 0L) Unfortunately, even though the gluino can decay dominantly to the bi, 
the latter (like the sl) does not decay to charged leptons and will not contribute to the 
like-sign dilepton signal. On the other hand, if there is a low-lying dn, one could still 
get a signal, although even this is suppressed by the fact that (3^ < 0.5. We do not 
go into this aspect any further since there is no compelling reason at the moment to 
postulate a light cIr. Thus we are forced to the conclusion that the like-sign dilepton 
signal is unlikely to be able to distinguish the stop solution to the HERA anomaly from 
the possible non-supersymmetric alternatives. However, a of mass upto 120-130 GeV 



is likely to be accessible to the third-generation leptoquark search at the Tevatron |2T 
(in the bb + missing Ex channnel) in the event of a luminosity upgrade. It might also be 
possible to access a part of the parameter space through cascade decays |22|,|2"3 



To summarize, then, we have investigated the possibility of observing like-sign dilep- 
tons at the Fermilab Tevatron in view of the possibility — no longer so remote, in view 
of the HERA anomaly — that there exists a ci squark of mass around 200 GeV. If there 
is also a gluino in the 200-350 GeV range, it might be possible to see such signals at 
the Tevatron, if not now, then in the near future, with the projected luminosity upgrade. 
Even with the present data sample it should be possible to explore a hitherto untouched 
part of the parameter space, at least for the scenario in question. Unfortunately, the alter- 
native SUSY solution of the HERA anomaly with a 200 GeV ti cannot be addressed well 
through this particular signal. Even so, much can be achieved by searching for like-sign 
dileptons in the present data sample and we would urge upon our experimental colleagues 
the importance of carrying out such searches with the highest priority. 



The authors would like to thank Manuel Drees and N.K. Mondal (DO Collaboration) 
for discussions. 
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